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Abstract 

An effective Hamiltonian for tlie localized spins in the one-dimensional Kondo 
lattice model is derived via a unitary transformation involving a bosonization 



'^ ' of delocalized conduction electrons. The effective Hamiltonian is shown to 

Q ' reproduce all the features of the model as identified in various numerical 

simulations, and provides much new information on the ferro-to-paramagnetic 

\f^ • phase transition and the paramagnetic phase. 
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The Kondo lattice model (KLM) describes the interaction between a conduction band 
and a half-filled narrow /-band, and is thought to capture the essential physics of some of 
the rare earth compounds ||1|. Although intensively studied for two decades, the KLM is 
still far from being completely understood. Even in the simple one-dimensional (ID) model, 
and with the conduction band less than half-filled, there are only two limits in which the 
behavior has been analyzed successfully; in the limit of vanishing conduction electron (CE) 
density, and for antiferro magnetic Kondo couplings J > 0, the /-electrons (/-spins) form 
a ferromagnetic (FM) ground-state 0; in the strong-coupling limit J -^ oo, and for any 
filling of the conduction band, the unpaired /-spins are again found to be FM [Q. The 
intermediate- to weak-coupling regime, away from half-filling but at finite CE density, has 
proved particularly difficult to analyze 0. 

From the known limiting behavior 0,^, together with a consensus of recent numerical 
simulations using the density-matrix renormalization-group, exact numerical diagonaliza- 
tion, and quantum Monte Carlo [^|-|§], a successful theory of the less than half-filled ID 
KLM will account for the following ground-state behavior of the /-spins: (i) At strong- 
to intermediate-coupling the unpaired /-spins are FM at all fillings and show behavior in 
accord with the strong-coupling expansion [^ . (ii) As the coupling is lowered, and for finite 
CE density, the system undergoes a transition to a paramagnetic (PM) state, with a filling 
dependent critical coupling in the weak to intermediate range, (iii) At weak-coupling, the 
system is characterized by a strong peak in the /-spin structure factor at 2kp of the CEs. 

In this Letter we derive an effective Hamiltonian i^eff from the ID KLM which repro- 
duces all the observed behavior in the intermediate- to weak-coupling regime, ifefr treats 
the /-spins exactly while the CEs are treated using bosonization techniques. The essential 
new ingredient in our work is an emphasis on describing delocalized CEs, as these are re- 
sponsible for the observed magnetic behavior of the /-spins. The problem of accessing the 
intermediate- to weak-coupling regime nonperturbatively is solved using a unitary transfor- 
mation. The effective Hamiltonian maps to the quantum random transverse-field Ising spin 
chain near the FM-PM boundary, and using extensive work on this interesting model by 



Fisher 0, we can obtain a vast amount of information on the transition and the properties 
of the model near it, as well as information on the PM phase. 
The Hamiltonian of the ID KLM is given by 

H = -tE(4s+i. + H.c.) + JES/,-Sc, (1) 

JO- j 

where t > is the CE hopping, S/j = lJ2a,a' c\j„(Ta^a'Cfj„,, Scj = | Ea,a' c]^cr<^,<^/Cj.^, and 
(T are the Pauli spin matrices. Fermi operators Cj^,Cj„ with subscript / refer to localized 
/-spins, those without refer to the CEs. We consider antiferromagnetic Kondo couplings 
J > and assume the conduction band filling n = Nc/2N < 1/2 with Nc the number of 
CEs and A^ the number of sites. 

From the strong-coupling expansion |^ , it is clear that the infinite J on-site spin-singlets, 
in which a CE is strictly localized with an /-spin, are magnetically inert: the strong-coupling 
FM only appears at large but finite J via CE hopping to neighboring unpaired sites, with 
a preferred spin orientation due to broken spin-singlet symmetry. The interaction identified 
in the strong-coupling expansion is the Zener double-exchange mechanism. This motivates 
us to introduce a delocalization length a > a {a the lattice spacing) which limits the mini- 
mum spatial spread of the CEs. The delocalization length models the qualitative difference 
between large J and infinite J behaviors, and has its physical basis in the energy gain for 
CE hopping to unpaired /-spins whenever t > 0. It relates to the average spatial spread of 
the CEs engaged in the double-exchange process. For example, the delocalization length in 
the one CE KLM corresponds to the effective spread of the spin polaron 0. For simplicity, 
a will be taken as an average applying uniformly to the CEs. It is important to empha- 
size that a limits only the minimum spread of the CEs and does not significantly affect the 
weak-coupling behavior, although it is essential in order to describe the strong-coupling FM. 

It is well-known that ID electrons may be represented using bosonization techniques. 
The Bose description is usually based on the Luttinger model due to its formal rigor, but 
this is not essential. In the present case it is essential not to use the Luttinger model, as 
will become clear. Two facts, peculiar to ID, form the basis of bosonization for realistic 



ID systems. The first is Tomonaga's observation [Q that the number fluctuation operators 
satisfy Bose-like commutation relations [pra{k), Pr'a'{k')] = Sr^r'^k-k'^ay rkL/2TT on a weak- 
couphng long-wavelength subspace, where the right-moving (r = +) and left-moving (r = — ) 
number fluctuations 

0<rp<7r/a 

with L = Na. The second is the fact that these number fluctuations generate the ID state 
space 10 . The main result from bosonization needed here is the representation of the Fermi 
site operators Cja in terms of the bosonic number fluctuations Praik). It is convenient to 
decompose the site operators into right- and left-moving components Cjo- = Y.r^^rju'- 

1 
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E e^^'^cfc. 



kp-±.<rk<kp+^ 



with kp = vrn/a, and where the momentum cutoff comes from Fourier analysis. In the 
Luttinger model the Bose representation may be formulated as an operator identity [jlO 



For the realistic system we must be satisfled with an approximate representation, but one 
which generates asymptotically exact results [|lT|. (The existence of the representation is 



guaranteed by the completeness of the Bose states.) In the thermodynamic limit, 

where the Bose flelds for u = p,a are deflned by 'ipv{i) = i{r[ / L)Y,k^o^^^'^°'[^+{k) ± 
i'^{k)]K{k)/k, with + corresponding to the number flelds ipy = (pu and — to the cur- 
rent flelds ipy = 6y. The charge and spin number fluctuations prik) = J^aPraik), and 
ar{k) = Yju '^ Pra{k) ■ Eq. (2) has the same form as in the Luttinger model but with one 
crucial difference. The even cutoff function K{k)^ satisfying A(A;) ^ 1 for \k\ < 1/a and 
A(fc) K, otherwise, is needed in the Bose flelds to ensure that delocalized CEs are de- 
scribed. The normalization factor N'{a) depends on both the cutoff and the cutoff function, 
and can only be determined asymptotically. Eq. (0) will of course fail if it is used to calculate 
number operators n^ja = clj^c^j^. In this case a Fourier expansion gives 
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dAMj) + ^Mj) + ^[Mj) + re^j)]} (3) 



to an additive constant. The separate form for the number operators is manifest also in the 
Luttinger model and is accounted for there with a carefully constructed normal ordering 
convention and a prescription for the correct taking of limits [|10 . 



To derive an effective interaction between the /-spins from the bosonized Hamiltonian 
(obtained by substituting Eqs. (0) and (0) in Eq. (|ID), it is sufficient to change to a basis 
of states in which the CEs are coupled to the /-spins. This is achieved using a unitary 
transformation with U = i{aJ /2'KVp) Yjj S'i, 0a{j), and where Vp = 2at sin{7m) . A variant 
of this transformation was first used by Emery and Kivelson for the single- impurity Kondo 
problem, and later generalized to the ID KLM [|I2|. The usage here is different; indeed 
the FM J^ term (see Eq. (|^) below), which U was designed to generate, is entirely absent 
in the previous work. The reason is that a Luttinger model bosonization will miss any /- 
spin effective interaction which is due to the non-local character of the CEs. Formally, in 
the Luttinger model the Bose fields (pu{j) and Huij) = —d^Oyi^j) are canonically conjugate 
and their commutator strictly vanishes unless they are at the same site. In our system 
the fields are smeared over a range a and their commutator is finite over roughly 27ra: 
[(f)y{j),Iiyi{^)] = 2i5yyjj{a) where Jj{a) = j^ cos{kja)Is?{k)dk. As examples, a Gaussian 
K{k) = exp(— a;^A;^/2) gives Jj{a) = (0r/2a) exp[— (ja/2a)^], and the Luttinger cutoff 
exp(— a;|A;|/2) gives J^j{a) = al\o? + (ja)^]. The Luttinger model (5-function is obtained 
by taking a — > in the last. The effect of this difference on the transformed Hamiltonian 
H = e~^ He^ is dramatic. Keeping all terms, 

+ JM\a) Y.ie'^^-^^'^'^'^S^ + H.c.]{cos[ir(j) - 0.(j)] + cos[2fcpja + 0,(j)]} 

i 
+ 2JM\a) E sin[ir(j) - 0.(j)] sin[2A;^ja + 0p(j)]5|, (4) 

where K{i) = —i{aJ/27iVF)J2j'['PuU)^^!^U')]^fj'- -^ condition for the derivation of Eq. (^ 
is that the cutoff be not too soft. 



The new term in Eq. (^) is the second. Since Si, is not transformed under U, it is 
immediate that the system is FM at intermediate-couphng at all fillings. The physical basis 
for the interaction is quite simple. A CE spread over more than one lattice site will carry 
the same spin over these sites. Due to the term J J^ji^rj^ — rirji) S^j in Eq. (|I]), this will tend 
to align the relevant /-spins. This interpretation also makes it clear that the interaction 
J^j{<y) is short-range provided a is finite. We may therefore approximate the FM term by its 
nearest-neighbor form —JcsJ2j ^fj^fj+iJ ^cs = (a^J^/27r^fir)j7'i(a). Although formally this 
term will give FM at strong-coupling as well, it is important to recall that the bosonization 
describes delocalized CEs. If J is too large then there will be significant CE localization and 
our approximation is less satisfactory. Note that it is in principle possible to include these 
effects as well with a less crude measure of CE delocalization and with the sum over j in 
the FM term restricted to sites containing unpaired /-spins only. Such alterations will not 
affect our conclusions, except to further support them. 

An effective Hamiltonian for the /-spins is obtained from Eq. (^ by replacing the CE 
Bose fields by their expectation values in the noninteracting ground-state. This step may be 
justified for the Bose charge-number field 4>p{j) by noting that at weak-coupling, which is the 
only regime where any of the fields affect Eq. (^), the charge structure factor is free electron 
like ip. For the spin fields there is less justification, though note that at weak-coupling 
these fields will be relatively smooth and will enter Eq. (^ as simple parameters. Thus 
while this approximation may affect the quantitative predictions of the theory, it would not 
be expected to affect the qualitative behavior. (Further evidence for this view was recently 
provided in a numerical simulation in which the same general behavior for the /-spins was 
seen with t — J interacting CEs [l^l-) The effective Hamiltonian is then 

j 

+ 2JAf\a) J2[cos K{j) + cos{2kF3a)]S'fj 
j 

+ 2JM'^{a) Y, sin K^j) sm{2kF3a)S}j (5) 

j 

and the spin directions have been reversed for later convenience. Eq. (Bl) is our main result. 



The remainder of this Letter is concerned with a brief analysis of ifgff to show that it gives 



all the required behavior. Details will be presented in a paper to follow ||14|| . 

To describe the destruction of the FM phase, the Sij, in K{j) may be replaced by 
their eigenvalues. K{j) is then a long-range object which counts the total Sf to the left 
of j and subtracts from that the total S*? to the right. (The effects of the non-Luttinger 
bosonization are not important here; [(f)u{j'),Gu{0)] -^ ^7rsgn(j') at large j'.) Near the 
FM phase boundary, and in the thermodynamic limit, it follows that K{j) ^ and any 
transition is described by the first two terms in Hes, with the second term responsible for 
spin-flips. For incommensurate fillings, cos{2kFJa) oscillates unsystematically with respect 
to the lattice. The large values cos{2kpja) ~ 1 which are responsible for spin flips, are 



then widely separated. Following analogous treatments in spin-glasses |jT^, this behavior 



is well-described by taking cos{2kpja) as a random variable. The factor multiplying S"?- 
in Eq. (^) is then replaced by hj, where hj is drawn independently from the displaced 
cos distribution p{h)dh where p{h) = (1/C7r){l - [{h/C) - 1]2}-V2 and C = 2JJ\f^{a). 
Note that fluctuations in the Bose charge-number fields 0p(j) offer further support for this 
interpretation. The behavior of the /-spins at and near the destruction of the FM phase 
is then governed by the quantum random transverse-field Ising spin Hamiltonian Hcnt = 
—JesY^j Sfj^fj+i ~'Yj ^j'^fj- Using extensive real space renormalization-group work on this 
model by Fisher (to whom we refer the reader for details), we determine the location of 
the quantum critical line describing the order- disorder transition at 

- = -——- sm{7m) . (6) 

t aji[a) 

The numerical predictive powers of Hes are restricted by lack of knowledge of A/'(a) . We 
would like to emphasize that such problems beset any bosonization description in which 
physical quantities are found to depend on this factor, and are not due to our particular 
bosonization. Accordingly, the coefficient of sin(7™) in Eq. (|^) is used as a fitting parameter 
to numerically obtained critical points 0-§|. A good fit is obtained with J/t = 2.5sin(7rn) 
as shown in Fig. 1. Note that this ignores any functional dependence of a on J or n. 



For the following discussion, it is convenient to introduce a measure of deviation from crit- 
icality 6 oc ln[2'n-^N"^{a)vp/a'^JJ'i{a)] [^, which for the obtained fit is 5 oc ln[2.5tsin(7m)/J]. 

The behavior described by i/crit is simply understood in terms of clusters of ordered 
/-spins. Reducing J from intermediate values in the FM phase, the infinite cluster charac- 
terizing strong FM is broken up into several large clusters as the quantum fluctuations hj, 
controlled by the spin-flip interactions, become stronger. The individual clusters are the spin 
polarons. The system is weakly ordered, and exists for —0.7 < 6 < with the boundary de- 
termined by Jcs = max{/ij}, as shown in Fig. 1. This is not a true transition line, but rather 



marks the onset of a Griffiths phase [|I6| characterized by singularities in the free energy over 
the whole range of 6. For small 6 the correlation length is ^ ~ (5~^, beyond which the system 
is ordered. The spontaneous magnetization Mq oc \S\^ with /? = (3 — -\/5)/2 ~ 0.38, while 
for small applied fields H the magnetization M{H) oc Mq[1 + 0{H'^^^^S\nH)]; the suscep- 
tibility is infinite with a continuously variable exponent. The mean correlation function 



{SfjSfjj^J — Mq oc |(5p'^(,^/a;)^/^e"''/^exp[— 3(7rx/,^)^/^] for x > ,^ and where the averaging 
is over p{h) []^. 

Further lowering J, we reach the true phase transition Eq. (j^). The correlation length is 
infinite, the magnetization M{H) oc | Inif |~^ for small H, and the mean correlation function 



Immediately below the critical line {S > 0), the system presents a weakly disordered 
Griffiths phase. The remaining clusters occupy a small fraction of the system length but 
"think" that they are still in the ordered phase; their magnetization 6^ per unit length is 
identical to Mq of the weakly ordered phase. These remaining rare clusters dominate the low- 
energy physics. The magnetization M(i^) oc 5>^{H'^^[S\n{l/ H) + const.] + 0[H'^^S\n{l/ H)]}; 
thus M{H) has a power law singularity with a continuously variable exponent 26; as in 
the weakly ordered phase the susceptibility is infinite. The mean correlation function 



decays less rapidly than in the ordered phase, but takes the same form {SfjSjj^^j.) oc 
(52^(^/x)5/6e-^/«exp[-3/2(7rx/0^/^] for x > ^ = 1/6"^ . According to i^crit, the weakly 
disordered Griffiths phase extends down to J = 0. However, as the disorder increases, the 
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third term in Hes is no longer negligible. At very low J, the last two terms in Hes will 
dominate; this corresponds to free spins in a field with dominant correlations at 2kF of the 
conduction band, and is responsible for the observed peak in the /-spin structure factor 
1^-0. No clusters remain. This strongly disordered conventional PM phase is indicated 
schematically in Fig. 1. 

In summary we have derived an effective Hamiltonian for the /-spins in the ID KLM 
which reproduces all the behavior seen in numerical simulations in the intermediate- to weak- 
coupling regime: (i) Hcs presents a FM phase at intermediate-coupling due to "forward" 
scattering by delocalized CEs, and is consistent with known limiting behavior PIQ. (ii) 
As J is lowered this phase is gradually disordered due to spin-flip interactions between the 
CEs and the /-spins. A sharp quantum order- disorder transition occurs to a PM phase 
at a critical coupling given in Eq. (j^). (iii) The backscattering interactions leave a residue 
correlation at 2k p in the /-spins at weak-coupling. 

This work was supported by the Australian Research Council. 
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FIGURES 

Fig. 1. Ground-state phase diagram of the ID KLM. The sohd (critical) hne is from Eq. 
(^) with 4:Ti'^N"^{a)/aJi{a) used as a fitting parameter to numerically determined points: 
square is density-matrix- renormalization-group data of a 75 site chain from Ref. [Q; dia- 
mond is the quantum Monte Carlo data for a 24 site system from Ref. |0; open circles are 
the exact numerical diagonalization data for the 8 site chain from Ref. |^. The dashed 
lines separate conventional strongly ordered (FM) /disordered (PM) phases from their weak 
(Griffiths phase) counterparts. 
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